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Heavy  Ion  Radiation  Effects  in  VLSI 


1.  INTRODUCTION 

Nonuniform  distribution  effects  show  up  dramatically  in  dosimetry  when  the 
dimensions  of  the  region  of  interest  are  approximate  to  the  ranges  of  the  electrons 
involved.  Decent  development  of  interface  dose  enhancement  phenomena  attests  to 
that  fact.  In  the  continuing  size  reduction  associated  with  microelectronic  circuitry 
a dimension  is  being  approached  in  which  uniform  dose  testing  will  be  inappropriate 
to  specify  response  on  hardness.  The  electrons  generated  by  heavy  ion  (for  example, 
cosmic  ray!  traversal  are  a case  in  point.  The  ion-electron  collisions  generate  a 
flux  with  ~1/K"  spectrum  and  maximum  energy  of  about  HI  keV,  The  range  in 
silicon  of  a 10  keV  electron  is  1.4  microns.  Thus,  when  microelectronic  devices 
reach  this  size,  local  dose  anomalies  can  occur.  Very  I.arge  Scale  Integration 
(VLSI)  crosses  the  threshold  into  this  domain. 

2.  DKVICE  SIZES  IN  VLSI 

l.Sl  technology  has  increased  the  density  of  circuit  elements  to  about  70,000 
per  chip  with  a device  size  of  about  1-mil  square.  Present  masking  techniques 
(optical)  lead  to  edge  acuities  of  2500  L,  This  means  that  the  intradevice  dimen- 
sions are  greater  than  a few  microns  (25  p - 1 mil).  As  a result,  the  energy  depo- 
sition profile  normal  to  an  energetic  heavy  ion  track  occupies  an  area  small  in 
(Deceived  for  publication  18  May  1978) 
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comparison  with  device  element  cross  sectional  areas.  Hence,  a single  ion  track 

cannot  influence  u device  of  l .SI  dimensions  except  in  the  most  contrived  and  extreme 

circumstances.  1 2 3 4 This  is  due  to  the  fact  that  the  energy  loss  mode  of  a heavy  ion 

is  electron  collisions.  That  collision  process  produces  a spectrum  of  electrons 

which  has  a few  members  of  energy  > 10  keV  and  the  range  of  a 10  keV  electron  in 
•> 

silicon  is  l.  4 microns.* 

In  consequence  of  the  foregoing,  a dramatic  effect  of  single  heavy  ion  energy 
deposition  will  not  become  apparent  until  device  sixes  approach  1 or  2 microns. 

This  situation  will  occur  with  the  advent  of  V 1 .SI.  The  projections  are  for  edge 
acuities  of  50  t that  could  lead  to  element  sizes  of  several  hundred  A and  devices 
about  l pon  edge. The  transition  is  made  possible  by  use  of  soft  x-ray  or  electron 
beam  masking  techniques. 


».  Sl/K.  AN1)  DOSK  D1STR1IUTIONS  NKAR  HKAVY  ION  TRACKS 


When  a fast  heavy  ion  penetrates  condensed  matter  some  of  the  orbital  electrons 
are  stripped  from  the  ion  and  a beam  of  such  ions  then  has  a mean  ionic  charge, 
t *■  /.  has  been  measured  widely  as  a function  of  energy  and  velocity  of  the  ion 
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and  is  given  by  Harkas 


./.a 


(1) 


The  ion  loses  energy  through  electron  collisions  and  generates  an  electron  flux  at 
a rate  given  by 


dn  2»  N x^  e*  1 , j 2 . , . . # ff  z , a)  , 1 , 

® -fir  V 


<21 


where 


dn 


» the  number  of  electrons  of  energy  u), 
u>  + die  generated  per  cm  of  ion  track  in 
matter  containing  N electrons/cm  , 


1.  Binder.  D.  et  al  (197M  1KKK  Trans.  Nuc.  Sci.  NS-22(No.  H:2675. 

2.  Ashley.  J.  C.  et  al  (1966)  HADC-TH-76- 125. 

3.  Holton.  W.  C.  (1977)  The  large  scale  integration  of  microelectronic  circuits. 

Set.  Am. 

4.  Harkas,  W.  H.  (1963)  Nuclear  Research  Emulsions.  Vol.  1 Academic  Press. 
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2 2 2 

u>ni  = 2 m c /3  y e,  m are  electron  charge,  mass. 

Integration  of  this  equation  yields  the  number  of  electrons  between  two  energies 


Wl*  W2* 


/ dn/du)  = 2nN(e  /me  ) 
W, 


2,2  m c 


2 , , , 2 
z If) 


(see  Reference  5). 

Thus,  the  integral  spectrum  falls  off  like  r[r.  Evaluation  of  the  integral  above 

w 

5 keV  yields 


/ dn/du)  = 2nN 
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{A) 
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which  for  Fe  cosmic  rays  with  (J  = 0.  5 in  silicon  gives  4.  3X  10  electrons/cen- 
timeter. 

Thus  in  a slab  of  Si  1500  A thick  only  6 electrons  are  liberated  with  energies 
greater  than  5 keV.  The  range  of  a 5 keV  electron  in  Si  is  4400  A so  the  electrons 
are  confined  to  a radius  < 4400  A (ignoring  diffusion). 

The  energy  loss  rate  in  the  target  is  obtained  from  the  Livingston- Bethe 
formalism 


aw  4.2 

, 4te  Z .. 
dE/dx  = r, N Z 


. 2m  v”  , . , ,,  2.  j 2 | 

Cn  — i C n(  1 - f)  ) - i) 


For  Fe'1^  cosmic  rays  this  can  be  ~30  Cle  V/ cm  so  that  in  a slab  1500  A thick 
~450  keV  could  be  dumped. 

A much  less  probable  but  more  dramatic  energy  loss  occurs  in  the  nuclear 
reactions  that  can  result  (p,n;  p,  2n;  p,e»,  etc.  ).  When  protons  or  a particles  are 
the  projectile,  these  reactions  lead  to  "star"  reactions  that  can  deposit  much  more 
ionizing  energy  locally  than  the  primary  proton. 

The  above  electron  generation  mechanisms  can  be  coupled  with  electron  trans- 
port models  to  produce  a relatively  detailed  picture  of  the  dose  and  electron  flux 
near  a heavy  ion  track.  Katz6  has  made  extensive  calculations  of  the  energy 
deposited  near  an  ion  track  with  the  following  result: 

5.  Mott,  N.  (1929)  Proc.  Hoy.  Soc.  124-425. 

6.  Kobetich,  E. , and  Katz.  R.  (1968)  Phvs.  Rev.  170,  No.  2. 
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l>  « K *“/«'“ 


that  is  the  energy  deposition  in  ergs /g  tit  radius  r from  a cylindrical  track  falls  off 
3 3 

as  l/r  and  is  proportional  to  t . where  from  above,  t is  the  mean  ionic  charge 

of  the  ion  or  cosmic  ray.  Both  t and  k are  functions  of  ft . Recent  measurements 

at  Hrookhaven'  have  verified  the  l/r"*  dependence  over  a range  of  1 to  300  nm 

(10  i to  3000  Al. 

2 

The  l/r  character  can  be  developed  approximately  as  follows; 


(11  The  initial  electron  spectrum  is  created  with  1/K“  probability 
0(K)  - 00/KJ. 

0.  7 

(3)  The  stopping  power  dE/dx  - k/E  ‘ at  low  energies, 

(3)  The  range  of  an  electron  at  low  energies.  R K K1'  ' . 


(41  The  electron  attenuation  rate  is  £(rl  - O^e  1 **. 


Therefore  the  total  energy  loss  rate  at  a radius  r is  given  by 
r 

max  ... 

dE«  - / *(r.  dE  dr  ■ 

0 

‘"max  «>„  . ,R  k/K0-7 

«.  - / _U_  e r/K  k/fc  dE  dr  . 

0 E 

Kmax  k^0  .r/kKi.7 

* - / — rr^r  e db  ur  . 


The  dose  (energy  deposited  in  a shell  dr  thick  divided  by  the  mass  of  the  shell) 

fc*  i «»  1 • 7 

. max  ;r/KK 

STf53r  * liTp  % { 0 {')  ' ■j.X’T <iK ' 

L *^Hmax  , *kK»0 

■ “3  e L " Tip"  • 

r 

2 

hence  since  is  of  the  order  of  1 M.  l)(r)  will  follow  as  l/r  for  the  first 

fraction  of  a micron. 


7.  Varma,  M.  N.  et  al  (1977)  Radial  dose  for  O*  in  Na«  Rad.  Res.  70:M1. 
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In  particular  the  dose  profile  in  quartz  for  Fe 
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cosmic  rays  (see  Figure  1) 
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F’igure  1.  Radial  Dose 
Distribution  Around  Heavy 
Ion  Track  Normal  to  Plane 
of  an  MOS  Device.  VLSI 
dimensions  are  approximate 


Thus,  the  picture  that  describes  the  generated  electron  energy  deposition  is  an 
intense  central  core  ~®k  in  radius  with  doses  ~107  rad  for  heavy  ions  and  10*’ 
rad  for  protons  and  mesons,  that  dose  falling  off  as  r “ with  a 100  rad  radius  of 
about  1 M in  silicon. 
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4.  STATISTICS  OK  THE  ENERGY  DEPOSITION  PROCESS 

When  a heavy  ion  traverses  a target  the  total  energy  loss  is  the  sum  of  the 
individual  electron  collision  losses.  If  the  number  of  collisions  is  large  and  the 
total  energy  loss  is  large  compared  to  the  maximum  loss  to  an  electron  in  a single 
encounter,  then  the  energy  loss  distribution  is  a gaussian.  If  not,  if  the  total 
collision  loss  is  comparable  to  the  maximum  single  loss,  then  the  distribution  is 
skewed  by  the  weight  of  improbable  but  large  energy  loss  single  collisions  and  a 
distribution  described  by  Landau,  Vavilov  or  Symon  results,  1 hus,  in  computing 
the  energy  loss  of  a heavy  energetic  ion  traversing  a thin  target,  one  must  take  care 
to  include  the  production  of  energetic  delta  rays  anti  not  simply  use  the  continuous 
slowing  down  approximation. 

Since  it  is  the  high  energy  electrons  that  are  rarely  produced,  one  must  use  the 
appropriate  statistics  to  describe  them.  The  generation  rate  of  electrons  of  energy- 
greater  than  some  E^  is  given  by  Eq.  (3).  If  one  takes  this  as  the  average  rate 
(delta  rays/ cm)  is 

P(n)  = i — — Mn  where  M is  the  Mott"'  result, 
nl 

For  example,  in  the  case  described  in  Section  2,  M = 6 electrons  (>  5 keV/  1500  A) 

P(6)  = 0.  13  for  example,  one  gets  6 electrons  only  13%  of  the  cases, 

P(  10)  = 0.  05  for  example,  one  gets  10  electrons  in  5%  of  the  cases. 

The  point  here  is  that  the  actual  yield  is  subject  to  considerable  variation  and 
using  only  the  mean  value  misrepresents  the  situation. 


5.  SOURCES 

Sources  of  fast  heavy  ions  are  the  naturally  occurring  fluxes  in  space  and  the 
man-made  accelerated  variety.  The  cosmic,  solar,  and  trapped  earth  protons  and 
ions  have  been  studied  and  measured  but  not  in  elaborate  detail.  The  proton  comp- 
nent  varies  due  to  magnetospheric  and  solar  activity  and  is  hard  to  describe  with 
a single  profile.  The  heavier  components  are  reasonably  stable  and  quite  energetic. 

For  radiation  effects  purposes  the  fluxes  and  energy  profiles  seem  well  enough 
known  to  make  calculations  on  the  number  of  events  to  be  expected  in  a satellite  due 
to  natural  background  of  cosmic  rays.  At  least  one  attempt  is  in  the  IEEE  literature. 

8.  Apparao,  M.  M.  V.  (1975)  Components  of  cosmic  radiation.  Topics  in  Astro 
Physics,  Vol.  11. 


The  flaxes  of  particles  can  be  made  into  5 groups:  protons;  a particles  (lie);  C.O,  N; 

Mg,  Si;  Fe,  Ni,  Co.  Protons  constitute  ~90%  of  the  total  flux;  He  about  9%.  The  LET 

o 

rates  range  up  to  30  GeV  cm  /g.  The  energies  of  the  most  energetic  10%  are 
greater  than  1 GeV /nucleon,  sufficient  to  penetrate  one  to  several  inches  of  alumi- 
num. 

Manmade  fluxes  can  be  generated  by  a variety  of  accelerators.  Proton  linacs 

capable  of  intense  beams  of  1 GeV  are  in  operation  and  the  Berkeley  Bevalac  has 

56 

recently  achieved  1 GeV/nucleon  Fe'  beams. 

The  cosmic  fluxes  can  be  evaluated  by  their  kerma  as  follows: 


.,  , , dN  dE  . . 

Hose  = k / ^ ^ dE 


in  particles/ cm^  MeV 


dE 


in  MeV  cm  / gm 


K = 1.  6 X 108  rad/g/MeV  . 

g 

This  has  been  done  for  tissue  with  the  following  results: 


Proton 

Z = 

1 

Dose  rate  rad/year 

4.  6 

He  ion 

Z = 

2 

3.  5 

M ion 

6 s 

Z s 

9 

1.  9 

LH  ion 

10  s 

Z s 

14 

1.3 

V H ion 

26  s 

Z s 

28 

1.  3 

EFFECTS 

The  effects  of  fast  heavy  ion  energy  loss  are  the  same  as  routinely  encountered 
in  radiation  hardening;  pseudo  permanent  damage  (for  example,  displacement,  oxide, 
and  insulator  charge  up)  and  transient  (induced  conductivity,  electron-hole  pair 
creation,  diffusion,  and  collection).  In  that  sense  the  description  is  unchanged. 

What  is  different  has  already  been  alluded  to  in  the  determination  of  the  sizes  of  the 
affected  areas.  The  following  differences  occur: 

(a)  Higher  doses  than  any  proposed  for  hardness  are 
achieved  in  the  central  core  of  an  ion  track.  Higher 
by  three  orders  of  magnitude  for  heavy  ions  (Fe), 

(b)  In  the  natural  background  case,  the  doses  are  delivered 
randomly  in  time  and  location  in  the  satellite, 

9.  Curtis,  S.  B.  (19)6)  IEEE  Trans.  Nuc,  Sci.  NS- 23- 1355- 1360. 
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(c)  Average  doses  due  to  run mtc  background  are  trivially  am. ill 
yrt  rarli  individual  event  which  would  cause  no  effort  in  SSI, 

May  trigger  1 SI  circuit  response  (see  Hinder  iffl  or  potential 
annihilation  of  a circuit  element  in  \ I .SI. 

(ill  The  dei  reuse  In  si *e  of  a circuit  element  leads  to  smaller 
p in  liability  of  heavy  Ion  lilt,  hut  the  purpose  of  \ I, SI  is 
to  Increase  the  total  number  of  circuits  m the  same  proportion, 
hence  the  "target"  area  remains  fixed, 

<«•>  Threshold  effects  (charge  required  to  flip-flop,  death  of 
biological  cell,  charge  release  needed  for  dielectric 
breakdown)  can  appear  which  would  never  be  seen  in  total 
dose  scenario.  As  an  example,  if  a local  energy  dump  ii: 
is  reipiired  for  threshold  event  then  AK  Am  is  local  dose 
required.  If  that  dose  wen  say  10  rad,  one  could  irradiate 
with  low  energy  electrons  or  photons  right  ..p  to  10**  rad  and 
see  no  events.  However,  every  heavy  ion  which  passes 
through  Am  (for  small  Am)  can  create  events  even  though 
average  total  dose  is  near  zero.  In  this  way  total  dose 
hardening  concepts  fail  and  must  be  replaced  with  a 
mterodosc  distribution  concept, 

(f)  These  effects  occur  without  a nuclear  warfare  scenario, 

(g)  These  effects  could  become  dominant  in  a star  wars  scenario 
(ion  beam  weapons,  satellite  satellite  encounters,  flooded 
phase  space  injection  with  fast  heavy  ions). 

INTER  ACTION  CROSS-SECTION  ANI>  ENEHliY  DEPOSITION 

Without  performing  detailed  calculations  one  can  still  scope  the  cosmic  ion 
interaction  problem  from  the  available  literature,  figures  -*  to  -1  show  coni|H<sites 
of  the  particle  spectra  for  protons,  carbon,  and  iron  and  their  respective  I I T's  as 
a function  of  energy /nucleon.  These  results  are  taken  from  deferences  !>  and  10. 

figure  fi  shows  the  integral  over  all  cosmic  species  for  I TIT  ' I I : T v s I I T. 
for  instance,  there  are  U50,  000  partieles/cm^’  year  with  l.K'T  ' 100  Me\  cm"  c. 
The  effect  of  ill)  g/em  of  shielding  (1.4  m.  All  is  also  shown. 


10.  Heinrich.  W.  (UUflHud,  Effects  .14: 14. 1 40. 


Ill 


Figure  -4.  I FT  and  Free  Space  Spectrum  For  Protons 


In  the  first  vertical  column  of  Figure  5 is  the  number  of  hits  on  assumed  VI. SI 
sensitive  volumes  in  a typical  chip,  0.2  cm  on  edge  containing  It'1  devices,  l'he 
area  chip  is  4 v 10  p (4  p“  per  device)  with  a depletion  depth  of  2 p.  The  average 
chord  length  for  a convex  device  (the  gate  region)  IS  2 \ 0.  2 pis  0.  3 microns, 
rhus.  the  mean  energy  deposition  by  an  ion  is  given  h\ 

- 4 ,3 

l.KT  \ 0.3  s 10  cm  \ 2.33  g/cm  for  silicon. 

I’hese  energy  dumps  are  shown  in  the  second  and  third  columns  of  Figure  5. 

The  computation  proceeded  as  follows: 

dF  2 3 — 

Me\  cm  g \ g/cm  V s cm  F Me)  . 


Also, 


,-tl 


— ergs  tOO  g • Me)  Mi  \ cm“/g  \ g/c’  V s cm  I*  rad  , 


s 4 V/A  for  convex  bodies  , 


l.  H v io‘H  t . vV  t •'  • 


114) 
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Figure  5.  Integral  LET  Spectrum.  Columns  show  hits  and 
energy  deposition  in  VLSI  sensitive  volumes  of  assumed  size 
(see  text) 

Thus,  the  average  dose  from  each  particle  to  any  convex  object  can  be  obtained 
from: 

1.  6 X 10~8  • (MeV  cm2/g)-  i = D rad  . 


A is  the  total  surface  area,  not  projected  cross  section.  For  the  example  above 

A = 2(2X  1 + 1 X 0.  2 + 2X  0.  2 )(i2  = 5.  2 X 10‘8  cm2  . 

2 

and  for  a maximum  LET  particle  = 6.  5 X 103  MeV  cm 

1.6X  10"8  X 6.5X  103  X jr  = 8 X 103  rad  . 

5.  2X  10 
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However,  as  previously  mentioned  this  represents  the  average  of  a distribution 
9 2 

that  varies  from  ~10‘  rad  to  10  rad  over  1 p radius  from  the  ion  track. 

The  maximum  energy  transfer  in  a single  collision  between  ion  and  electron  is 


E = 4 <me/Mion>E.on 


For  Fe  at  1 GeV  this  is  35  keV.  Hence,  those  energy  dumps  * 70  keV  can 
be  expected  to  show  the  Vavilov  distribution. 


a SIMULATION  OF  HEAVY  ION  EFFECTS 

Experimental  simulation  of  the  effects  of  cosmic  ion  irradiation  can  be  achieved 
to  a large  extent  (not  displacement  effects  and  restricted  to  thin  targets,  t < few 
microns)  by  the  electron  beam  in  a scanning  electron  microscope.  Those  effects 
deriving  from  low  energy  electron  collisions  such  as  electron-hole  pair  production 
and  subsequent  trapping  would  be  reproduced  well.  Others,  some  more  lethal, 
would  not.  (For  instance  the  creation  of  a pin  hole  punch  through  leading  to  a gate- 
substrate  short.  ) In  Table  1 the  comparisons  are  indicated  for  various  important 
parameters  of  simulation. 

Figure  6 shows  the  differential  free  space  flux  of  large  LET  particles.  One 
can  use  the  graph  to  determine  the  relationship  between  SEM  current  and  the  cosmic 
ray  flux  at  each  LET  rate.  Figure  7 shows  the  radial  dose  profiles  from  a SEM 

operated  at  10  and  20  keV  as  computed  by  Chadsey ' ' and  the  previously  mentioned 

6 7 

work  of  Katz  and  of  Varnia. 

12 

Galloway  and  Poitman  have  indicated  the  pertinent  parameters  for  an  SEM 
as  a device  to  deliver  uniform  exposure  to  LS1/MOS  devices,  but  nonuniformity  is 
the  characteristic  of  radiation  effects  in  the  small. 


11.  Chadsey,  W.  L.  (1973)  NAD  Crane  Technical  Kcport  7024-C74-69. 

12.  Galloway,  K. , and  Poitman,  P.  (1977)  IEEE  Trans.  Nuc.  Sci. 

NS-24(No.  6): 2066. 
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Figure  7.  Comparison  of  Radial  Distributions  For  Heavy 
Ions,  SEM  Electron  Beams,  Theory,  and  Experiment 

9.  CONCLUSIONS 

Nonequilibrium  dosimetry  effects  become  apparent  when  regions  of  interest 
(devices)  have  dimensions  approximate  to  the  ranges  of  the  electrons  involved. 
Total  dose  concepts  fall  in  that  domain.  VLSI  radiation  effects  testing  and  analysis 
will  require  elevated  sensitivities  amongst  the  radiation  effects  community  to  cope 
with  these  phenomena.  The  cosmic  ray  interactions  in  space  applications  are  a 
case  in  point. 
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The  numbers  suggest  an  enhanced  vulnerability  to  damage  and  device  mode 
change  in  VLSI  size  devices.  The  S15M  simulation  offers  a partial  solution  to  the 
assessment  problem. 
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